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We present preliminary results on a study about the thermal variation of the charged and neutral 
pion masses to one loop, analyzing their electromagnetic difference, in the context of Chiral 



Perturbation Theory with two flavours, as well as using a light resonance model. We find that 
f^h the pion mass difference increases for, at least, low and intermediate temperatures, unlike the 

chiral limit decreasing result. The axial-vector mixing arising from chiral restoration smooths the 
Debye- screening temperature increase. Taking into account further corrections due to axial and 
vector resonances, dominated by a; and p particles respectively, does not change significantly the 
ChPT prediction. 
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1. Introduction 

SU(2)-Chiral perturbation theory (SU(2)-ChPT) is one of the most powerful techniques to 
describe the low energy sector of the strong interaction. It is an effective field theory based on the 
spontaneous breaking of the chiral group, SUv(2) x SUa{2) — > SUv(2), and provides a systematic 
and model-independent scheme for calculating low energy observables at zero [1, 2, 3] and also at 
finite temperature [4, 5], which is very important to compare with the results provided by lattice 
calculations [6, 7] below the pseudocritical temperature (180-200 MeV). 

The terms of the effective Chiral Lagrangian are constructed as an expansion in even powers of 
a meson energy scale (external momentum, meson mass) which is to be compared finally with the 
chiral scale, A ~ 1 GeV, that separates in an approximate way the region of high energy processes 
from the low energy ones. 

In order to implement electromagnetic (EM) effects, the vector subgroup SUv(2) must be 
explicitly broken by adding a photon field and a quark charge matrix via the external source method, 
giving rise to new terms in the Chiral Lagrangian [8, 9, 10] and introducing the need of a change 
in the chiral counting scheme in order to accommodate the electric charge, e, which is formally 
achieved considering e 2 = 0(p 2 /F 2 ), being F the pion decay constant in the chiral limit. 

We present here preliminary results for the leading order corrections to the real part of the 
self-energy of neutral and charged pions propagating in a thermal bath. Our analysis could be 
useful when comparing neutral and charged pion distributions observed in Heavy Ion Collisions. 
Actually, recent experimental data are able to measure and compare charged and pion distributions 
very accurately [11, 12], showing that their difference is negligible within errors at very low pj, 
where our approach would be applicable. Besides, there are finite temperature lattice analysis 
measuring pion screening masses [13] which should not differ much from the pole masses analyzed 
here below but not close to the transition. EM interactions are customarily not included in lattice 
simulations and therefore an study of their effect in different thermal quantities is interesting in 
order to support or not such approximation. Finally, the analysis has also formal interest, since the 
pion mass difference is connected to the axial-vector spectral function difference at T = [14] and 
at T / in the chiral limit [15]. Since non-zero quark masses invalidate by definition the use of the 
soft pion limit leading to such connection, the analysis for non-zero masses and finite temperature 
is pertinent and the role of resonances can be established by the use of effective models. The 
detailed analysis, results and discussion will be the subject of a forthcoming work [16]. 

The expression for the pole of the pion propagator establishes the dispersion relation as cus- 
tomary, as Pq= p 2 + M 2 + L K (po , p) where p is the three-momentum of the external particle, M n is 
the tree level mass of the pion, and £^ is the &(p A ) pion self-energy, which encodes all the interac- 
tions with the medium and will depend separately on po and p as a result of the Lorentz Symmetry 
breaking due to the choice of the reference system associated with the thermal bath. If EM isospin 
breaking is considered, the dispersion relation is different for charged and neutral pions already at 
tree level. The real part of L x gives information about the dispersion relation and the pion velocity, 
whereas the imaginary part is related to the absorption in the medium. Since we are interested only 
in calculating masses, we will focus our study in the real part and will make use of the position of 
the pole of the in-medium pion propagator, i.e, our modified masses of charged and neutral pions 
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Figure 1: 1-PI diagrams contributing to the self-energy of a charged pion in SU(2)-ChPT to leading order. 
Diagrams for neutral pions are the same removing those in which photons are present. 

will be the solutions of 

Ml i (p)-Ml i -ReZ Ki (p 2 =p 2 +Ml h p) = 0, i = {0,±}, (1.1) 

where the self-energy is calculated perturbatively in ChPT and therefore is on-shell, with a residual 
thermal three-momentum dependence. We will define, as customary, the pion masses in the static 
limit p = 0. 

2. Pion masses in SU(2) ChPT 

The effective chiral lagrangian up to fourth order in p (a meson mass, momentum, temperature 
or derivative) including electromagnetic interactions proportional to e 2 is given schematically by 
Jz?e/y = J£ p 2 +e 2 + Jz? p 4 +e 2 p 2 +( ,4. The second order lagrangian corresponds to the familiar non-linear 
sigma model plus the addition of the gauge coupling of mesons to the photon field via the covariant 
derivative, and an additional term proportional to a low-energy constant C compatible with the 
e 7^ symmetries of the QCD lagrangian [8, 17]: 

F 2 

^ p i +e i = — tr[D^D^U + 2B ^ (U + U f )] +Ctr[QUQU*] , (2.1) 

where U(x) = exp[£<J>/F] G SU(2), being <I> the pion matrix field; and the covariant derivative is 
defined through Dp = dp + iAp [Q, •] with Ap the photon field. ^# and Q are the quark mass and 
charge matrices defined by Jt = mid and Q = (e/3)diag(2, -1), and C = {M 2 ± -M^F 2 jle 2 . The 
fourth order lagrangian consists of all possible terms compatible with the QCD symmetries to that 
order, including the EM ones, and can be found for SU(2)-ChPT, for example, in [9]. We will not 
have to be worried by the set of EM and non EM low energy constants introduced by ^ p 4 +e 2 p 2 +e 4, 
for we will be mainly interested in thermal contributions and the fourth order lagrangian is just used 
to absorb the T = divergences coming from the loops. The complete zero temperature one-loop 
values for the pion masses including EM corrections can be found for instance in [9, 18]. 

The self-energy of a charged pion to leading order consists of all the 1 -particle irreducible 
(1-PI) diagrams modifying the pion propagator at 0(p 4 ) in SU(2)-ChPT , shown in Fig.l. The 
neutral pion does not feel the electromagnetic field to this order, so the corrections to its mass do 
not involve photon fields (the rest of the diagrams are the same although the vertices are different). 

The photon-tadpole diagram Fig. 1(c) is proportional to the photon mass and therefore vanishes 
at zero temperature, while pion tadpoles Fig. 1(a) and the photon-exchange diagram Fig. 1(d) are 
finite and chiral scale-independent once regularized and combined with diagram Fig. 1(b). 
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The pion tadpoles, charged or not, give rise to g\ (M, T) thermal functions defined asg{(Mf,T) 

27t 2 



!-2 Jo dp^--p^ — , with Ep = p 2 +Mf and j3 = T 1 . For the neutral pion mass we get 



M 2 =M 2 n0 {T = 0) 



\+ Vl \gim ± j)-- gl {M^x, 



(2.2) 



The above correction due to the EM tree masses correction represents numerically only a slight 
deviation from the M 2 ± = M 2 limit, the neutral pion mass still increasing with temperature. 

The charged pion mass receives contributions from two different sources: pion tadpoles and 
virtual photon diagrams, which also can be of two different types: a photon-tadpole contribution 
(Fig. 1(c)) which grows like e 2 T 2 , just in the same manner as a Debye or screening mass of the 
electric field in a thermal bath [19]; and a more complicated structure of one-photon exchange 
(Fig. 1(d)) that depends on the external pion momentum. Putting together all the pieces we get 

M 2 a 4p) =M 2 n± -AZe 2 gl (M 2 n± ,T) + l -e 2 T 2 +ReL rEx (p), (2.3) 

being p the external pion three-momentum, Z = C/F A and LyEx(p) the self energy contribution 
from the one -photon exchange on the mass shell [18, 16]. The differences between considering 
the static limit (p = 0) and values of the three-momentum compatible with the temperature range 
considered are tiny [16], so that we will restrict ourselves to this limit to define the pion masses. 

Results for the charged and neutral masses separately, and for their difference can be seen in 
the left plot 1 of Fig. 2. We have used physical masses for the pions instead of the tree level masses 
for the numerical results, since the difference between these and the tree level ones is encoded in 
higher order corrections. Despite the different sign of the various terms, the screening contribution 
(quadratic in T) turns out to dominate the pion mass difference, which, as the masses separately, 
increases only slightly near the critical temperature with respect to its T = value. When T grows 
above the applicability range of these ChPT calculations the mass difference decreases, as should 
be expected since expansions in M n ± /T — > should coincide with the chiral limit calculation of 
previous works that found a T-decreasing behaviour for the difference [15]. For physical masses 
and realistic temperatures, our result is qualitatively different from the chiral limit one. 

Actually, and following this last line of action, we may wonder what would our results become 
when taking the limit of temperatures (i) much greater than the masses and the external momenta 
(which means that we have to set the masses inside the loops to zero), and (ii) sizeable to the 
momenta, k, running inside the loops, i.e T ~ k 3> m,p; where p here is the external momentum of 
the pion. Expanding our expressions leads to the result given in [15], M 2 ± —M 2 = M 2 ± ^1 — ^ J + 

\e 2 T 2 , which serves us as another consistency check. 

As commented above, our low temperature analysis allows to work with a slightly different 
chiral limit, in the sense that we can still assume m = but considering e ^ 0, even inside the 
loops. In the right plot of Fig. 2 we show our calculation both in this latter limit, and also considering 
m j£ 0, e 0; to be compared with those appearing in the chiral limit result [15]. The screening-like 
terms, always increasing with T and inherent to the thermal bath, compensate for the terms coming 



'We point out that Fig.4 and Fig.5 shown in [18] do not correspond to the analytical expressions for the masses of 
the charged and the neutral pion given there, which are correct. 
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Figure 2: Left: Charged (red,dashed line) and neutral (blue, solid line) pion masses in the static limit to 
leading order in SU(2)-ChPT for non zero tree level pion masses. Right: Different results for the charged- 
neutral pion mass difference, (a, solid line) corresponds to our results in the chiral limit keeping corrections 
e ^ for the tree level charged pion mass inside the loops, (b, dot-dashed line) corresponds to the same 
calculation with m ^ and e =^ also inside the loops; and the full dashed line is the result given in [15]. 



from the sum rule relating the axial and vector spectral functions [14], which should decrease 
applying chiral symmetry restoration arguments [21]. The sum rule is not valid for non-zero quark 
masses, but one can still identify the contributions from the vector and axial spectral functions from 
a resonance exchange model, as we will briefly describe in the next section. 

3. Resonance contributions 

The results for the pion mass difference in ChPT is quantitatively reliable only in a temperature 
range up to, at most, 100-150 MeV. One way of increasing the predictive quality of the calcula- 
tions and check the validity range of the ChPT result -at the risk of losing the model-independent 
character of the theory- is to use a model [17, 20] that implements the couplings with the vector res- 
onance J PC = 1 and the axial / = 1 ++ . This allows to calculate the corrections to the charged 
pion mass (neutral pions do not couple to this interactions in the model) allowing to extend the 
calculations up to intermediate energies (~1 GeV) [22]. 

The couplings with the spin-1 resonance fields relevant for the leading order calculation of 
the charged pion self-energy are diagrammatically expressed in Fig. 3. It is convenient to note that 
we still have the contributions coming from charged and neutral pion tadpoles and virtual photon 
diagrams. Since resonance fields only couple to charged pions, the neutral pion mass at leading 
order in this resonance model is the same as that obtained in the previous section. 

As commented before, at T = and in the chiral and soft pion limit, it was shown [14] that 
the pion mass difference can be parametrized exactly in terms of the difference of the vector-axial 
spectral functions, so the corrections to this formula must be proportional to M n , as was shown in 
[23] using ChPT at zero temperature. The modifications to the expression of Das et al. for non-zero 
quark masses and temperatures are discussed in [16]. 

Although there are more realistic calculations making a non-narrow width treatment of the 
resonances [22], we will suppose that the spectral functions pv{s) and Pa{s) are of the form 
Pva{s) = Fy A 8(s — My A ), where Fy, Fa, My and Ma stand for the vector (axial) decay constants; 
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Figure 3: 1-PI diagrams contributing at leading order to the charged pion-self energy, p and ai particles 
are represented by double and dashed lines, respectively. The relevant vertices including charged pions, 
resonances and photons are drawn as grey boxes. 
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Figure 4: Charged and neutral pion masses in the static limit to leading order in the light resonance model 
(left plot), and their difference (right plot, solid line) compared to the leading order ChPT calculation (right 
plot, squares). 

and vector (axial) masses, respectively. We will use the numerical values for this constants given 
in [22]. The thermal modification of these spectral functions amounts to leading order to the pion 
tadpole contributions and accounts for the p — a \ mixing in the thermal bath [21]. Other corrections 
-such as the modification of the p width showing up in the dilepton spectrum- are suppressed in 
the pion mass difference for the temperature range considered here. 

The results for the charged pion and the neutral masses in the static limit, as well as the pion 
mass difference in this model compared with that calculated in SU(2)-ChPT to leading order, are 
shown in Fig.4. The first we observe is that contributions coming from the interaction with reso- 
nances are only activated at temperatures much higher than those at where low energy effects start 
to be important (which correspond to the ChPT results), as we would have expected since, up to 
T ~ 200 MeV, resonances are suppressed. The sub-leading resonance corrections are typically of 
&(T 1 M\lM 1 p ai ). Due to this very same fact, the contribution corresponding to the exchange of an 
ai and a photon is much smaller than those in which the p interacts (form factor and tadpole-like 
diagrams). It turns out that the inclusion of resonant contributions does not change significantly the 
ChPT result (extrapolating the thermal range in which results are realistic, we get approximately 
5% of discrepancy at 250 MeV). Of course, the lack of an appropriate coupling constant for per- 
turbative calculations implies the fact that there may exist other higher order diagrams which can 
be numerically relevant at finite temperature, and we must not exclude these resonances to have a 
more important role even in this intermediate energy region, as they have at zero temperature [14]. 
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